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SUMMARY 


This  is  the  Phase  III  of  AF  Contract  F  33615-79-C-5143:  Leadtime 
Variability  in  Inventory  Requirements  Projections  —  Final  Technical  Report 
and  Sumnary.  The  report  summarizes  the  first  two  technical  reports  and  the 
results  of  the  cost  benefit  analysis.  The  study  consisted  of  three  phases, 
and  this  is  the  third  phase.  The  purpose  of  the  study  was  to  find  out  whether 
leadtime  variability  had  a  significant  influence  on  stock  safety  levels. 

We  found  that,  indeed,  it  has. 

Currently,  AFLC  revises  its  leadtime  according  to  the  following  model: 

"Lt  =  Lt-i  • 

A 

where  Lfc  is  the  estimate  of  the  required  leadtime  during  period  t  and  Lt-1 
is  the  actual  leadtime  observed  during  the  preceding  period.  This  model  is 
inappropriate  because  it  implicitly  assumes  that  leadtime  for  the  next-period- 
ahead  is  deterministic,  whereas  it  is  stochastic.  Cohen,  in  a  similar  study 
for  the  U.  S.  Army,  makes  like  observations  (Cohen,  1979). 

The  first  phase  of  the  study  was  concerned  with  fitting  distributions 
to  leadtime  data  of  high-intensity  items.  The  results  showed  that  procurement 
leadtime  can  be  fitted  very  well  by  the  lognormal,  the  gamma,  the  normal,  or 
the  Weibull.  The  gaimia  distribution  has  very  often  been  postulated  for  the 
distribution  of  leadtime  (Burgin,  1972,  for  example)  and  our  results  support 
this  assumption. 

The  second  phase  was  to  be  a  simulation  study  and  an  analysis  to  examine 
the  impact  on  inventory  control  levels  of  the  fitted  distributions  of 
leadtime.  The  simulation  model  was  to  have  been  provided  by  AFLC,  but  in  a 
February  5,  1980,  meeting  at  WPAFB  with  Major  Paul  Gross  of  AFBRMC  and  Gloria 


Picciano,  Oiann  Lawson,  and  Carol  Hawks  of  LORRA,  it  was  agreed  that  it  would 
be  more  beneficial  to  do  independent  simulations  as  well  as  look  into  several 
important  issues  not  specified  in  the  original  contract.  The  oral  agreement 
was  embodied  in  a  memo  of  understanding  received  on  March  7,  1980.  This  memo 
is  provided  as  Appendix  A  to  this  report. 

Addressing  the  several  issues  raised  by  the  memo  of  understanding,  we 
first  conjecture  that  the  "skewed  normal"  distribution  reported  upon  by  Demmy 
(1979,  p.  IV-13)  in  his  analysis  of  forecast  error  is  likely  produced  by 
faulty  analysis.  We  support  this  conjecture  by  simulation  and  by  an  analysis 
the  reader  will  find  in  section  3.  A  good  forecasting  methodology  should 
produce  normally  distributed  forecast  errors. 

When  leadtime  and  demand  are  stochastic,  leadtime  demand  would  be  a 
compound  distribution  or  a  convolution.  We  knew  the  distribution  of  leadtime 
for  the  items  in  the  sample,  but  we  did  not  have  the  corresponding  distribu¬ 
tions  for  demand.  Consequently,  we  assumed  that  daily  demand  at  a  depot 
is  Poisson-distributed.  We  then  obtained  via  simulation  the  compound 
distribution  of  leadtime  demand  for  a  few  representative  items.  And  we 
succeeded  in  fitting  the  right  tails  of  these  simulated  leadtime  distribu¬ 
tions  to  the  Laplace.  This  result  is  important  because  it  verifies  an  important 
assumption  in  the  work  of  Presutti  and  Trepp  (1970).  With  Laplace  leadtime 
demand,  Presutti  and  Trepp  have  already  worked  out  the  optimization  models. 

We  also  found  that  the  normal  distribution  is  a  good  approximation  to  the 
distribution  of  leadtime,  but  that  the  Laplace  is  better.  Optimization 
techniques  concerning  the  normal  have  been  extensively  worked  out  (see  Brown, 
1967).  Section  5  of  this  report  deals  at  some  length  with  the  distribution 
of  leadtime  demand. 


Section  6  presents  results  of  a  cost-benefit  analysis.  Using  a  scenario 
at  WR-ALC  with  shortage  factor  X  *  660,  requisition  size,  R  =  4,  and  essentiality 
parameter,  Zi  =  0.5,  we  simulate  the  current  system  versus  the  Presutti-Trepp 
model  IV  (1970,  pp.  248-249)  optimal  system.  We  find  that  the  current  system 
is  wanting,  the  present  average  service  level  being  79.6%  versus  an  optimal 
level  of  86.9%.  To  bring  the  system  to  optimality  for  our  sample  requires  a 
one-time  investment  of  about  $750,000  in  safety  stocks  and  an  increase  of 
about  $250,000  in  annual  operating  costs.  With  an  annual  procurement  of  one 
billion  dollars,  our  sample  represents  0.5%  of  that  procurement.  Thus,  if  the 
sample  is  representative,  AFLC  would  require  a  one-time  investment  of  $150 
million  in  safety  stocks  and  an  increase  in  annual  operating  cost  of  $50  million. 

It  is  also  evident  that  the  present  mix  of  stock  levels  is  inappropriate, 
and  that  the  policy  of  maintaining  an  aggregate  safety  level  stock  not 
exceeding  a  two-month  supply  is  not  sound.  The  reader  will  discern  these 
details  in  section  6. 
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1 .  DESCRIPTION  OF  THE  STUDY 

In  recent  years,  the  Air  Force  Logistics  Command  (AFLC),  with  an 
inventory  of  economic  order  quantity  (EQQ)  of  over  two  billion  dollars,  has 
encountered  inventory  support  problem  apparently  because  of  leadtime  varia¬ 
bility.  The  purpose  of  this  study  is  to  determine  whether  leadtime  varia¬ 
bility  has  a  significant  impact  on  inventory  support  planning  and  control  at 
AFLC. 

Using  Air  Force  leadtime  data  on  high-intensity  items,  we  performed  the 
runs  test  for  randomness  (Siegel,  pp.  52-56).  The  leadtime  data  were  random 
except  for  a  few  instances  of  suspected  "pencil -whipping."  With  random 
leadtime  data,  we  could  then  fit  statistical  distributions  to  them.  This 
was  accomplished  by  means  of  the  Kolmogorov-Smirnov  (KS)  goodness-of-fit 
test  (Conover,  1971,  pp.  293-298),  fitting  the  data  to  the  exponential, 
the  ganma,  the  normal,  the  Weibull,  and  the  lognormal.  The  results  indicated 
that,  in  general,  the  lognormal,  the  gamma,  the  normal,  and  the  Weibull 
could  be  postulated  as  good  fits  to  the  leadtime  data. 

The  impact  of  leadtime  variability  on  stock  levels  can  be  known  exactly 
if  we  identify  the  distribution  of  leadtime  demand.  To  do  that  we  need  to 
know  not  only  the  distribution  of  leadtime  but  also  of  period  demand.  Since 
we  had  no  demand  data  except  item  monthly  demand  rates,  we  assumed  that  daily 
demand  is  Poisson-distributed.  Then  with  the  best  fit  for  leadtime  and 
with  demand  Poisson,  we  simulated  the  distribution  of  leadtime  demand  for 
some  representative  items  in  the  sample.  The  results  indicated  that  the 
right  half  of  the  simulated  distribution  can  be  fitted  very  well  to  the 
Laplace.  This  is  quite  an  important  finding  since  AFLC  would  like  to  use 
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the  Presutti-Trepp  Model  IV  (1970,  pp.  248-249)  which  assumes  that  leadtime 
demand  is  Laplace-distributed. 

The  final  stage  of  the  study  was  to  conduct  a  cost-benefit  analysis  of 
current  and  proposed  stock  levels.  This  was  accomplished  with  the  aid  of 
data  received  from  LORRA  on  May  23,  1980.  With  that  data  and  using  Changes 
2  and  3  of  AFLC  Regulations  57-6  (Department  of  the  Air  Force,  29  December 
1978  and  22  June  1979),  we  were  able  to  describe  the  current  system.  And 
by  using  the  Presutti-Trepp  Model  IV,  we  were  also  able  to  generate  a  proposed 
optimal  system. 
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2.  SUMMARY  OF  THE  FIRST  TWO  TECHNICAL  REPORTS 

According  to  the  contract,  the  first  technical  report  (Phase  I)  was  to 
be  delivered  on  December  31,  1979,  the  second  (Phase  II)  on  March  31,  1980. 
The  Phase  I  report  was  to  address  the  statistical  distributions  of  leadtime, 
the  Phase  II  report  the  AFLC  simulation  model.  But  because  leadtime  data 
were  obtained  in  two  subsets  of  sizes  n  *  16  and  n  =  46,  the  first  in  mid- 
October  1979  and  the  second  toward  the  end  of  November  1979,  and  because  of 
the  memo  of  understanding  (see  Appendix  A),  the  first  two  reports  under- 
standaDly  did  not  conform  precisely  to  the  letter  of  the  contract. 

The  Phase  I  report  (27  pages)  advanced  the  concept  of  leadtime  demand 
as  a  compound  distribution  of  period  demand  and  leadtime.  It  also  reported 
the  results  of  goodness-of-fit  tests  on  the  first  subset  of  sixteen  items. 
Before  doing  the  goodness-of-fit,  the  data  were  verified  to  be  random. 

The  Phase  II  report  (96  pages)  completed  the  fitting  of  statistical 
distributions  to  the  leadtime  data,  addressing  the  second  subset  of  forty-six 
items.  The  Phase  II  report  was  made  of  seven  technical  appendices  in  order 
to  facilitate  future  research  and  to  serve  as  tutorials  on  the  statistical 
distributions  considered.  Perhaps  a  very  important  Appendix  in  the  Phase  II 
report  is  the  one  that  theoretically  examines  the  influence  of  leadtime 
variability  on  buffer  stocks  and  service  levels,  with  examples  provided 
on  selected  items. 
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3.  THE  DISTRIBUTION  OF  FORECAST  ERRORS 

In  his  analysis  of  forecast  errors,  Denmy  (1979)  arrived  at  distributions 
that  are  very  much  skewed  to  the  right,  distributions  called  "skewed  normal" 
by  LORRA.  The  purpose  of  this  section  is  to  show  that  when  there  is  a  mix 
dominated  by  declining  demands  for  the  many  and  sharply  increasing  demands 
for  the  few,  then  the  method  used  by  Demmy  is  very  likely  to  produce  dis¬ 
tributions  of  forecast  errors  that  are  skewed  to  the  right.  We  go  on  to 
validate  our  assertion  by  means  of  a  set  of  experiments.  But  first  we  must 
say  that  data  must  be  stationary  and  also  random,  that  is  uncorrelated, 
before  distributions  are  fitted  to  them.  This  was  evidently  not  in  the 
case  with  Demmy ‘s  forecast  errors.  We  concentrate  on  his  Table  IV-2  (Demmy, 
1979,  p.  IV-13).  The  distributions  he  produced  were  based  on  approximately 
22,500  items.  For  each  of  these  items,  he  calculated  a  simple  average 
and  the  MAD  of  the  demands  for  the  first  eight  quarters,  that  is,  FY  71-72 
(Demmy,  1979,  pp.  III-l  to  1 1 1-3  and  1 1 1-9  to  III-ll).  The  simple  average 
became  the  forecast  for  each  quarter  in  FY  73-75.  Thus  the  forecast  errors 
of  quarters  9,  10,  .  .  .  ,20  were  the  actual  demands  lor  the  quarter  in 
question  minus  the  simple  average  for  the  first  eight  quarters.  Demmy 
then  divided  this  forecast  error  by  the  historical  MAD  for  that  item  in 
order  to  obtain  standardized  scores  and  in  order  to  be  able  to  aggregate 
the  22,500  odd  forecast  errors  for  each  period.  He  called  these  standardized 
scores  Z9,  Z10,  ....  Z20. 

If  the  demands  for  these  items  were  stationary,  then  the  distribution 
of  forecast  errors  would  have  been  symmetric  if  not  normal.  But  these 
demands  were  not  stationary,  and  it  seems  that  most  were  declining  with 
time.  Consequently,  the  majority  of  the  forecast  errors  were  negative. 
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which  explains  the  negative  bias  in  forecast  errors  that  Demmy  obtained 
(Demmy,  1979,  p.  III-16).  With  declining  demands  it  is  obvious  that  the 
absolute  value  of  the  bias  would  increase  with  time:  this  is  precisely 
what  Oemmy  reports. 

In  summary  we  conjecture  that  tne  "skewed  normal"  distributions  obtained 
by  Demmy  are  merely  a  result  of  wrong  statistical  analysis.  Any  forecasting 
method  worth  its  salt  should  yield  normally  distributed  forecast  errors. 

Had  Demny  accounted  for  the  trends  in  demand,  his  forecast  errors  would 
have  been  normally  distributed;  and  to  validate  our  assertion,  we  perform 
a  set  of  experiments. 

A  Set  of  Experiments 

Eighty  items  are  considered  in  the  set  of  experiments.  The  demand 
function  of  every  item  can  be  described  by  the  following  expression: 

0 ( t )  =  D(t-l)  +  C  +  X  ,  (3.1) 

where  0(t)  =  demand  in  time  period  t, 

C  =  a  constant, 

and  X  =  a  normally  distributed  random  variable  with  zero 
mean  and  std.  dev.  equal  to  ox- 

Each  item  belongs  to  one  of  the  four  classes  described  below: 

Class  I  :  C  is  negative. 

Class  II  :  C  is  zero. 

Class  III  :  C  is  positive. 

Class  IV  :  C  is  positive,  but  has  a  value  of  C'  during  the  first 

eight  time  periods,  and  a  value  of  C"  during  the  next 
twelve  time  periods.  In  addition  C"  is  greater  than  C‘. 

Each  item  has  a  unique  combination  of  C  and  a*.  For  each  item  twenty  demand 

values  were  generated  by  a  mechanism  consistent  with  (3.1). 


Class  I  items  show  declining  demand  with  a  linear  trend.  Class  II 
items  exhibit  no  trend,  whereas  Class  III  items  show  increasing  demand 
again  with  a  linear  trend.  Class  IV  items  also  exhibit  increasing  demand 
but  the  trend  changes  abruptly  after  the  eighth  period. 


Forecasting  Technique  1:  Trend  not  removed.  The  average  demand  in  the 
first  eight  periods  was  the  forecast  for  each  of  the  next  twelve  periods. 
The  standardized  forecast  error  Z(t)  was  computed  as  follows: 


where 

and 


Z(t) 


F(t) 


0(t)  -  F(t) 
iMAO 


forecast  for  period  t 


t  -  9,  10,  ...,  20, 


8 

2  D(t) 

JM _ 
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MAD  *  mean  absolute  deviation  of  demand  during  the  first  eight 
periods. 


Forecasting  Technique  2:  Trend  removed.  The  first  difference  S(t)  =  D(t)  - 
D(t-l),  was  computed  for  the  first  eight  periods.  Values  of  0(0)  were  generated 
by  the  appropriate  mechanism.  The  average  value  of  S(t)  in  the  first  eight 
periods  was  a  forecast  of  S(t)  in  each  of  the  next  twelve  periods.  Let  the 
forecast  of  S(t)  be  FS(t),  and  the  forecast  of  demand  be  SF(t). 


Then  FS(t)  *  ±  S^q-  —  -  f  t  =  9,  10, 

and  SF(t)  =  D(t-l)  +  FS(t)  ,  t  =  9,  10 . 20 


20; 


The  standardized  error  SZ(t)  is  given  by 


Experiments  and  Results  —  Sample:  Type  A.  Type  A  samples  consisted  of 
items  from  Classes  I,  II,  and  III.  Each  sample  had  eighty  items.  When  the 

proportion  of  Class  I  items  in  the  sample  was  higher  than  that  of  Class  III, 

the  distribution  of  Z(t)  was  found  to  be  skewed  to  the  right  with  a  negative 

mean.  For  higher  values  of  t,  the  mean  shifted  further  away  from  zero.  The 

distribution  of  SZ(t),  however,  was  symmetric  with  a  mean  very  nearly  equal 
to  zero. 

As  more  and  more  Class  III  items  were  substituted  for  Class  I  items, 
the  distribution  of  Z ( t }  became  increasingly  symmetric.  When  the  proportion 
of  Class  III  items  exceeded  that  of  Class  I  items,  the  distribution  of  Z(t) 
became  skewed  :o  the  left  with  a  positive  mean. 

The  distribution  of  SZ(t)  remained  symmetric  in  spite  of  changes  in  the 
mix  of  items. 

Experiments  and  Results  —  Sample:  Type  B.  Type  B  samples  consisted  of 
items  from  all  four  classes.  A  typical  sample  had  the  following  composition 
40  Class  I  items 

20  Class  II  items 

10  Class  III  items 

10  Class  IV  items. 

An  example  of  these  classes  of  trend  is  given  in  Figure  3.1. 

For  this  sample  the  Z(t)  distribution  was  skewed  to  the  right  with  a 
long  right  tail.  This  may  be  immediately  seen  in  Figure  3.2;  and  for 
comparison  we  reproduce  in  Figure  3.3  Demmy's  Z9  distribution.  The  tail 
vanished  when  Class  IV  items  were  dropped,  though  the  distribution  remained 
skewed  to  the  right.  The  SZ(t)  distribution  was  also  found  to  have  a  long 


Figure  3.1.  Tour  Classes  of  Trend  Mixed  Together 
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right  tail.  With  the  elimination  of  Class  IV  items,  the  tail  disappeared 
and  the  distribution  became  symnetric  around  a  mean  that  was  very  close  to 
zero,  as  can  be  seen  in  Figure  3.4. 


In  conclusion,  the  assertion  by  LORRA  that  the  distribution  of  forecast 
errors  is  highly  skewed  to  the  right  does  not  seem  valid.  It  is  contrary 
to  theory,  and  is  most  likely  based  on  undue  confidence  in  a  very  faulty 
forecast  method. 
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4.  ABC  ANALYSIS 

Because  of  the  large  number  of  items  in  the  inventory,  it  is  sometimes 
convenient  to  classify  them  according  to  the  annual  dollars  of  usage;  and 
even  though  we  are  dealing  only  with  high  intensity  items,  it  would  still 
be  useful  to  classify  these.  For  the  fifty-six  items  in  the  sample,  we  find 
that  the  total  annual  usage  is  $4,868,000,  comprising  about  155,000  units 
of  demand.  Of  these  9.78%  account  for  54.34%  of  the  annual  dollars.  We 
call  these  Class  A  items  because  they  should  receive  the  most  attention. 

We  similarly  see  that  25.64%  of  the  number  of  units  account  for  only  9.53% 
of  the  annual  dollars.  We  call  these  Class  C  items  because  managing  them 
can  be  kept  as  simple  as  possible.  Between  A  and  C,  we  have  Class  B  items. 
For  our  sample,  these  comprise  64.58%  of  the  number  of  units  but  account 
for  36.08%  of  the  annual  dollars.  These  results  are  taken  from  Table  4.1, 
and  although  arbitrary,  the  classification  is  consistent  with  the  ABC 
classification  scheme  advocated  in  the  literature.  See,  for  example.  Brown 
(1967,  pp.  23-24)  and  Peterson  and  Silver  (1979,  pp.  71-73).  The  ABC 
classification  is  also  consistent  with  the  Supply  Management  Grouping  Code 
(SMGC),  even  though  the  latter  does  not  appear  to  be  as  useful.  If  we 
inspect  Table  4.1,  for  example,  we  see  that  SMGC  code  M  accounts  for  79.46% 
of  the  annual  dollars  but  32.51%  of  the  annual  units;  SMGC  Code  P  accounts 
for  20.38%  of  the  annual  dollars  but  66.66%  of  the  annual  units;  and 
SMGC  Code  T  accounts  for  0.16%  of  the  annual  dollars  and  0.83%  of  the 
annual  units.  In  terms  of  the  span  of  control,  it  should  be  clear  that 
an  ABC  classification  is  superior  to  the  SMGC.  Furthermore,  it  would  be 
easy  to  do  tradeoffs  between  inventory  investment  and  service  level  using 
ABC  curves  as  Herron  (1976)  shows. 
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In  Figure  4.1  is  plotted  the  cumulative  percentage  of  items  along  the 
horizontal  axis  versus  the  cumulative  percentage  of  total  dollar  usage  along 
the  vertical  axis.  Figure  4.1  makes  it  more  convenient  to  reclassify,  if 
necessary. 

According  to  Brown  (1967,  p.  23),  the  distribution  of  annual  dollars 
of  demand  is  lognormal.  This  is  indeed  the  case  for  our  data.  The  Lilliefors 
(L)  test  of  lognormality  yields  a  Dmax  =  0.1097,  less  than  the  critical 
statistic,  Lq  gg  =  0.386//5F  =  0.1184.  Hence  we  accept  the  hypothesis  of 
lognormality  for  the  annual  dollars  of  demand.  The  fit  to  lognormality  is, 
in  fact,  quite  good.  (The  parameters  of  the  lognormal  are  u  =  10.369  and 
a  =  1.213.) 

In  a  similar  vein,  we  also  fit  the  number  of  annual  units  demanded  to 
the  lognormal  distribution.  This  yields  Dmax  *  0.0890,  as  compared  to  the 
critical  LQ  g5  =  0.1184.  Thus  we  also  accept  the  hypothesis  of  lognormality 
for  the  number  of  items  demanded.  (The  parameters  of  this  lognormal  are 
u  =  7.263  and  a  =  1.157.) 

The  1980  Data 

The  data  used  to  generate  the  Figure  4.1  were  obtained  from  L0RRA  on 
May  23,  1980.  In  inspecting  these  data,  it  became  clear  that  the  sample 
of  sixty-two  we  originally  used  has  now  been  reduced  to  one  of  fifty-six 
items.  This  was  because  five  of  the  items  were  no  longer  high-intensity, 
and  one  item  was  inadvertently  duplicated  in  the  original  sample. 

For  the  sake  of  future  research  and  in  order  that  our  results  may  be 
replicated,  we  organize  the  data  received  on  May  23,  1980  in  Tables  4.2  and 
4.3.  These  tables  account  for  fifty-six  items  only;  we  have  kept  the  original 
item  numbers  reflecting  the  listing  in  the  Phase  I  and  Phase  II  reports. 
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Figure  4.1.  ABC  Classification  of  the  Sample  Items--1980  Data 


Source:  Table  4.1 
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Table  4.3.  DATA  OF  MAY  23,  1930  BY  ITEM  NUMBER 
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Table  4.2  lists  the  items  by  ascending  order  of  stock  number;  Table  4.3 
lists  these  items  by  ascending  order  of  item  number.  Note  in  Table  4.3  that 
these  items  were  deleted:  item  4  (FSN  3021 01866582JH :  Latch), 

5  (FSN  1 56000572261 60H:  Rod  Assembly),  6  (FSN  5305001 117850AB:  Screw  Spec), 

14  (FSN  589500172114:  Antenna),  16  (FSN  7045008479020),  and 

49  (FSN  4935006506352AB:  Cable  Assys  -  2).  Item  16  is  in  fact  item 

50  (FSN  7045008479020);  the  other  items  are  no  longer  high-intensity. 

The  1977-73  Data 

We  had,  prior  to  obtaining  the  1980  data,  played  detective  with  1977-73 

order  quantities  and  unit  prices  originally  provided  with  leadtime  data.  We 

thus  generated  the  annual  units  demanded  as  well  as  the  dollar  values  of 

demand  for  1977-78.  For  the  sake  of  comparison  with  1980  data,  we  provide 

Table  4.4  and  Figure  4.2.  Table  4.4  is  the  analogue  of  Table  4.1,  and 

Figure  4.2  is  the  analogue  of  Figure  4.1.  In  Table  4.4,  we  see  that  5.63ri 

of  the  top  units  account  for  63.09"  of  the  dollars;  56. 8:*  of  the  intermediate 

units  account  for  34.79"  of  the  dollars;  and  37.52."  of  the  bottom  units 

account  for  only  2.12"  of  the  dollars.  In  1977-78,  the  sample  accounted 

for  an  annual  usage  of  S6, 940, 000  and  comprised  399,523  units  of  demand. 

Thus,  the  demand  for  these  items  has  fallen,  on  the  average,  between  1977 

and  1980.  An  item-by-item  comparison,  using  Tables  4.1  and  4.4,  shows  that 

a  majority  have  declined  in  demand,  but  that  for  a  few,  demand  has  sharply 

increased.  This  is  consistent  with  our  conjectures  in  Section  3. 

The  distribution  of  annual  dollars  of  demand  in  Table  4.4  is  also 

lognormal.  The  Lilli efors  (L)  test  of  lognormal ity  yields  a  D  of  0.07017. 

max 

This  is  much  less  than  the  critical  statistic  Lq  ^  *  0.886/ *57  *  0.1125. 
Hence  we  accept  the  hypothesis  of  lognormal ity,  the  parameters  of  the 
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lognormal  being  p  =  10.595  and  a  =  1.210.  Furthermore,  the  number  of  annual 
units  demanded  is  also  lognormal  with  parameters  p  =  7.717  and  a  ■  1.407. 

The  L  test  yields  a  0max  of  0.1097,  less  than  the  critical  Lg  g5  =  0.1125. 

It  should  be  observed  that  the  analysis  cf  the  1977-78  imputed  data 
is  based  on  the  original  sample  of  size  sixty-two  which  contains  one  redundant 
item  (#16)  and  five  items  (-4,  5,  6,  14,  and  49)  that  are  no  longer  high- 
intensity.  Since  item  16  plays  a  very  minor  role,  the  effect  of  its  presence 
is  negligible. 
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5.  THE  DISTRIBUTION  OF  LEADTIME  DEMAND 

We  may  think  of  leadtime  demand  as  a  random  sum  of  demands  that  are 
independently  and  identically  distributed.  Thus  leadtime  demand  may  be 
written  as 

X  =  D1  +  D2  +  ...  +  D.  +  ...  Dl,  i  =  1,  2,  ....  L  ,  (5.1) 

where  D  and  L  are  random  variables  denoting  demand  and  leadtime.  Thus  leadtime 
demand,  X,  may  be  thought  of  as  a  mixture,  while  leadtime,  L,  is  the  mixing 
distribution.  More  specifically,  f(X)  may  be  said  to  be  a  compound  distribution 
with  G( L)  being  the  compounding  distribution  (Ord,  1972,  pp.  64-66). 

It  can  be  shown  (Drake,  pp.  109-112)  that  for  the  structure  (5.1) 

E(X)  -  E(L)  .  E(D)  ,  (5.2) 

and 

V*(X)  =  E(l)  .  V(D)  +  [E(D)]2  .  V(l)  ,  (5.3) 

the  star  denoting  the  variance  of  leadtime  demand  with  variable  leadtime. 

If,  on  the  other  hand,  laadtime  is  constant  at  L,  then  E(X)  would  be 
as  in  (5.2)  but 

V ( X )  =  E(L)  .  V ( D )  .  (5.4) 

We  can  immediately  see  the  influence  on  safety  stocks  if  we  begin  to 
consider  the  variability  of  leadtime.  For  the  same  safety  factor,  k,  this 
increase  would  be  in  the  ratio 
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E(L)  V(0)  +  [E(D)r  V(L 
^  E (L) V(D 


A  1/2 


•AS}-  e(d>) 


,  where  VMR  denotes  the  (5.5) 
variance  to  mean  ratio 


,  for  VMR(D)  =  1  and  large  VMR(L)  . 


A  similar  analysis  appears  in  Appendix  G  of  the  Phase  II  report  (Hayya, 

March  31,  1980,  pp.  88-94).  An  explicit  treatment  of  the  impact  on  safety 
stocks  will  be  considered  in  the  next  section. 

Theory  of  Compound  Distributions 

It  would  be  easy  if  compound  distributions  were  readily  recognizable: 
sometimes  they  are,  sometimes  they  are  not.  Hadley  and  Whitin  (1963,  p.  117) 
have  shown  that  where  the  procurement  leadtime  is  gamma-distributed,  with 
parameters  a,  3,  and  if  a  Poisson  process  with  mean  Xt  generates  demands  with 
units  being  demanded  one  at  a  time,  then  the  distribution  of  leadtime  demand 
is  a  negative  binomial  with  parameters  a  +  1 ,  3/(3  +  X).  Burgin  (March  1972) 
has  treated  the  case  with  demand  normal  and  leadtime  gamma.  There  is  other 
work  (for  example:  Sherbrooke,  1968;  Qrd,  1972;  Bott,  1977),  but  the  leadtime 
distributions  are  too  complex  for  our  present  purposes. 

In  the  absence  of  data  on  daily  demands,  the  assumption  of  Poisson  demands 
may  be  appropriate.  Furthermore,  the  fitting  of  statistical  distributions  to 
leadtime  data  supports  the  notion  that  leadtime  is  gamma-distributed.  Con¬ 
sequently,  we  can  assume  that  leadtime  demand  is  a  negative  binomial.  Better 
yet,  we  can  use  simulation  in  order  to  see  whether  we  can  fit  the  distribution 
of  leadtime  demand  for  our  items  to  the  normal  or  the  Laplace.  Particularly, 
if  we  can  fit  the  distribution  of  leadtime  demand  to  the  Laplace,  we  can  take 
advantage  of  the  models  developed  for  the  Air  Force  by  Presutti  and  Trepp  (June  1970) 
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Simulation  Experiment  1:  Poisson  Demands;  Best  Fit 
for  Leadtime— 1977-78  Demands 

We  had  no  demand  distribution  data  at  the  writing  of  this  report;  and 

until  May  25,  1980,  we  had  no  1980  daily  demand  data.  Consequently  we  made 

the  assumption  that  daily  demand  at  the  depot  >s  Poisson-distributed,  and,  at 

the  beginning,  we  estimated  mean  daily  demand  by  averaging  the  quantities 

ordered  over  the  two  most  recent  years  of  data:  1977  and  1978.  Then  we 

estimated  the  annual  dollars  of  demand.  These  ranged  from  $6,000  for  item 

no.  3  (Integrated,  FSN:  5962004537739)  to  $2,825,000  for  item  no.  37 

(20  mm  GunBa,  FSN :  1005007879802).  On  the  other  hand,  the  standard  deviation 

★ 

of  leadtime  demand,  a.,  ranged  from  29.33  for  item  no.  12  (Antenna,  FSN: 

5821002694508)  to  15,385.47  for  item  no.  23  (PAD-ASSY  1,  FSN:  1095009120256). 

★ 

We  could  see  by  inspection  that  dollar  value  of  demand  and  o i  are  hardly 
related.  A  calculation  of  the  linear  correlation  coefficient  yields  r  *  0.256, 
barely  significant  at  the  0.05  level  (the  calculated  value  of  the  test  statistic 
is  t(60)  =  2.055).  It  may,  however,  be  possible  to  fit  an  exponential  function 
to  the  two  variables. 

★ 

A  cross-tabulation  of  annual  dollars  of  demand  versus  produces  the 
frequencies  given  in  Figure  5.1.  From  the  cross-tabulation  in  Figure  5.1, 
we  choose  five  items  for  a  simulation  experiment.  These  are: 


m  no. 

Noun:  FSN 

* 

£1 

Annual  $  *  s 

23 

Pad  Assy-1:  1095009120256 

Very  high 
(15,385) 

High 

($173,000) 

37 

20  mm.  Gun  Ba:  1005007879802 

High 

(7,423) 

Extremely  high 
($2,325,000) 

39 

Guide  Rolle:  1005007545293 

Low 

(384) 

Very  low 
($11,000) 

50 

Lead  Tape:  7045008479020 

Very  high 
(10,530) 

Very  low 
($16,000) 

53 

Adapter:  1 5600094987JH 

Moderate 

(1346) 

Moderate 

($79,000) 

Annual  $'s  (in  thousands) 


Lead time  Demand--! 977-78  Data 
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The  best  fit  for  the  leadtimes  and  corresponding  parameters  is  given  in 
Table  5.1.  Table  5.2  gives  the  mean  and  variances  of  leadtime  demand.  With 
Poisson  demands  and  with  leadtimes  as  described  in  Table  5.1,  we  generate  by 
simulation  distributions  of  leadtime  demand,  using  samples  of  size  n  =  30  and 
n  *  100.  We  then  fit  these  distributions  to  the  normal  and  to  the  Laplace. 

The  results  are  given  in  Table  5.3:  we  see  in  that  table  that  we  are  not 
successful  in  fitting  the  normal  or  the  Laplace  for  the  entire  distribution 
of  leadtime  demand;  we  are,  however,  quite  successful  in  fitting  the  normal 
and  the  Laplace  to  the  right  half  of  the  distribution  of  leadtime  demand, 
that  is,  for  k  >  0. 

Simulation  Experiment  2:  Poisson  Demands;  Best  fit 
for  Leadtime--l98U  uemands 

Since  demands  declined  in  general  from  1977  to  1980,  we  verified  the 
results  of  the  previous  simulation  using  1980  demands.  Using  samples  of  size 
n  *  30,  we  again  simulate  daily  Poisson  demands  and  the  leadtime  distributions 
given  in  Table  5.1. 

The  means  and  variances  of  leadtime  demand,  using  1980  data,  are  given 
in  Table  5.4;  it  is  the  analogue  of  Table  5.2.  The  goodness-of-fit  of  the 
1980  leadtime  data  (actually  1980  demands  and  1975-78  leadtimes)  to  the  normal 
and  the  Laplace  are  also  verified  as  seen  in  Table  5.5.  It  is  seen  in  Table  5.5 
that  the  fit  is  better  in  the  right  tail  of  the  distribution;  the  fit  improves 
dramatically  beyond  k  =  1.  This  is  what  we  are  looking  for. 
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TABLE  5.1 

BEST  FIT  FOR  LEADTIME 


Item  no.  FSN 

Best  Fit 

Parameters 

23:  1095009120256 
(Pad  Assy-1) 

Normal 

u 

=  384 . 57 ;  a 

=  272.00 

37:  1005007879802 
(20  mm.  Gun  8a) 

Gamma 

a 

=6.09;  3  = 

75.19 

39:  1005007545293 
(Guide  Rolle) 

Lognormal 

u 

=  4.90;  c  = 

0.41 

50:  7045008479020 
(Lead  tape) 

Gamma 

a 

=  5.63;  3  * 

26.60 

53:  1 560009492087JH 
(Adapter) 

Gamma 

a 

=  3.40;  3  » 

91.74 
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Reject  the  hypothesis  of  good  fit  at  the  ninety-five  percent  significance  level. 
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TABLE  5.4 


MEANS  AND 

VARIANCES  OF  LEADTIME  DEMAND 
1980  DATA 

Item  no. 

Expected 

Daily 

Demand  E(Di) 

Expected 
Leadtime 
in  Days 

E(Li ) 

Expected 

Leadtime 

Demand 

E(X) 

23 

37.4 

384.6 

14,384.04 

37 

24.2 

459.3 

11 ,115.06 

39 

3.2 

145.3 

464.96 

50 

2.9 

149.6 

433.84 

34 


Variance  and 
Standard  Deviation 
of  Leadtime  Demand 

mh  q(x) 

103,500,305.7;  10,173.51 
20,283,043.47;  4503.67 
27,486.32;  195.79 
33,881.76;  184.07 
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4.1 


312.1 


1,279.61  434,781.98;  659.38 
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6.  COST-BENEFIT  ANALYSIS 


In  doing  the  cost-benefit  analysis,  we  compare  the  current  system  at 
AFLC  with  an  optimal  system.  Because  of  arguments  presented  in  the  previous 
section,  we  choose  the  Presutti-Trepp  Model  IV  (1970,  pp.  248-249)  as  the 
optimal  system.  One  simulation,  using  1980  data  and  with  X  *  660,  Z^  s  0.5, 
ai  =  0.23,  is  presented  in  Table  6.1,  which  is  sorted  by  FSN  for  convenience. 
And  based  on  the  same  parameters  and  upon  an  understanding  of  the  present 
procedures  at  10RRA,  Table  6.2  presents  a  description  of  the  current  system. 
How  each  column  in  these  tables  is  calculated  may  be  respectively  seen  in 
Appendices  8  and  C  of  this  report. 

The  question  concerning  the  mix  of  items  at  an  ALC  and  how  that  mix 
would  change  is  answered  by  comparing  the  corresponding  columns  in  Tables  6.1 
and  6.2.  The  mix  will  change. 

It  may  be  worthwhile  to  compare  the  current  system  with  the  equivalent 
optimal  system.  This  comparison  is  presented  in  Table  6.3  in  terms  of 
several  attributes:  service  level,  number  of  backorders,  fill  ratio,  value 
of  safety  stock,  and  so  on. 

It  may  also  be  worthwhile  to  study  the  behavior  of  the  optimal  system 
for  different  shortage  factors,  \.,  and  different  item  essentiality, 

Table  6.4  presents  this  behavior  for  Z^  =  0.5  and  for  X  =  660,  600,  500,  400, 
300,  200.  Table  6.5  gives  a  similar  comparison  for  Zi  =  1 

From  Tables  6.4  and  6.5  we  produce  the  exchange  curves  in  Figures  6.1 
and  6.2.  These  exchange  curves  should  give  an  idea  of  the  additional  invest¬ 
ment  in  safety  stock  required  to  reduce  the  number  of  backorders  or  to  improve 
the  service  levels.  The  reader  should  keep  in  mind  that  these  tables  and 
exchange  curves  are  based  on  a  sample  of  size  n  =  56.  This  sample  represents 
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an  annual  procurement  of  $4.9  minion.  Since  total  annual  procurement  is 
one  billion  dollars,  our  sample  represents  about  one-half  of  one  percent 
of  total  procurement.  Consequently,  the  numbers  in  this  section  must  be 
multiplied  by  200  in  order  to  fathom  the  impact  on  the  total  inventory  system. 
Returning  to  Table  6.3,  for  example,  we  see  that  to  bring  the  system  to 
optimality  would  require  a  one-time  additional  procurement  for  safety  stocks 
of  $727,400  x  200  -  $150  million.  We  can  also  see  that  the  annual  operating 
cost  will  increase  by  $240,234  x  200  =  $50  million.  This  is  the  price  to  pay 
in  order  to  reduce  backorders  by  6d63  x  200  =  1,400,000  units  per  year. 

In  Figures  6.1  and  6.2,  we  see  that  the  current  policy  itself  is  not 
on  the  optimal  curve.  Thus,  AFLC  can  presently  reduce  the  number  of  back¬ 
orders  without  any  additional  investment  in  safety  stocks.  Conversely, 
with  the  current  number  of  backorders,  AFLC  can  reduce  its  investment  in 
safety  stocks.  We  illustrate  by  means  of  Figure  6.2.  At  present,  the  total 
number  of  annual  backorders  for  the  sample  is  about  10,000;  the  total  safety 
stock  is  about  S600.000.  With  the  same  investment  in  safety  stock,  we  can 
immediately  reduce  our  annual  backorders  to  about  5750  units  per  year.  This 
translates  to  a  reduction  of  backorders  for  the  total  system  of  about 
(10,000  -  5750)  X  200  =  850,000  units  per  year.  Thus  the  number  of  backorders 
can  be  almost  halved  by  incorporating  leadtime  variability  and  by  adhering  to 
the  Presutti-Trepp  optimal  model. 
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7.  LIMITATIONS,  ASSUMPTIONS,  RESEARCH  QUESTIONS, 
METHODS  AND  ANALYSES  USED 


Limitations 

The  uninvolved  reader  may  find  fault  with  a  study  of  this  kind,  but 
considering  the  limited  scope  originally  envisaged,  that  reader  may  perhaps 
discern  the  amount  of  work  that  was  nevertheless  put  in  it.  We  are  dealing 
here  with  inventories  whose  worth  runs  into  the  billions  and  whose  management 
is  necessarily  very  complex.  Consequently,  it  would  be  vain  to  claim  that 
this  modest  study  will  resolve  all  the  problems  of  shortages  encountered  by 
AFLC. 

Assumptions 

Since  we  did  not  have  daily  demand,  but  only  demand  rates,  it  was 
convenient  to  assume  that  daily  demand  was  Poisson-distributed.  LORRA  had 
supplied  us  with  daily  shipment  data  on  two  items,  but  shipment  data  are  not 
demand  data.  The  daily  shipment  data  do  not  appear  to  conform  to  the 
Poisson  distribution. . 

We  had  also  assumed  that  leadtimes  and  demand  rates  were  statistically 
independent.  Lacking  simultaneous  demand  and  leadtime  data,  it  was  impossible 
to  validate  that  assumption. 

Finally,  we  implicitly  assumed  that  the  sample  of  fifty-six  high-intensity 
items  is  representative  of  the  population  of  such  items.  In  fact,  it  was  not 
a  random  but  a  convenience  sample. 
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Research  Questions 

In  a  study  of  this  kind,  the  questions  raised  by  the  research  overwhelm 
the  answers  provided,  and  this  study  is  no  exception.  Also  the  questions 
raised  depend  very  much  on  the  reader's  point  of  view.  From  our  point  of 
view,  there  are  numerous  questions,  and  as  an  example  we  ask  only  two: 

1.  Why  the  incongruities  in  the  mathematics  of  the  current  system? 

In  calculating  the  safety  factors,  AFLC  relies  on  a  Presutti-Trepp 
formula  which  in  turn  is  based  on  leadtime  demand  being  Laplace- 
distributed.  But  in  calculating  the  quarterly  MAO,  AFLC  pretends 
that  leadtime  demand  is  normally  distributed.  One  effect  of  this 
incongruity  is  in  the  calculation  of  the  standard  deviation  of 
leadtime  demand,  a,  from  MAD.  We  see  that  for  the  normal, 

o  »  1.25  MAD  ,  (7.1) 

whereas  for  the  Laplace, 

3  =  MAD  .  (7.2) 

Thus,  with  Laplace  leadtime  demand,  we  should  use  (7.2)  rather  than  (7.1). 

2.  What  is  a  proper  definition  of  fill  ratio?  Does  fill  ratio  address 
the  percentage  of  requisitions  filled  or  the  annual  fraction  of 
demand  satisfied?  Using  the  latter  as  a  definition  of  theoretical 
fill  ratio  does  not  produce  results  that  are  compatible  with  actual 
fill  rates,  as  may  be  seen  in  Table  7.1.  Of  course,  the  calculation 
of  theoretical  fill  ratio  in  Table  7.1  does  not  consider  customer 
priority. 
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TABLE  7.1 

CURRENT  FILL  RATES:  REPORTED  VS.  THEORETICAL 


Fill  Rates 


Item  no. 

FSN 

Reported^ 

2 

Theoretical 

1  -  bt/q 

9 

5895004451349 

81 .90 

99.0% 

99.6 

8 

5895009497160 

63.06 

85.0 

93.8 

7 

5960004009106 

98.79 

74.5 ' 

74.9 

3 

5962004537739 

67.80 

90.4 

92.3 

10 

5895000894403 

98.75 

80.6 

91.7 

2 

5841 000738392 

100.00 

95.9 

98.74 

15 

589001167508 

73.29 

82.4 

93.33 

11 

5831008803563 

99.76 

69.0 

67.4 

13 

5821009906461 

100.00 

100.0 

100.0 

12 

5821002694508 

37.28 

91 .2 

98.2 

50 

7045008479020 

99.41 

93.6 

99.7 

51 

3040006211345 

17.09 

93.3 

98.2 

52 

1670007970253 

72.57 

100.0 

100.0 

53 

1560009492087 

95.08 

81.3 

87.6 

28 

1095008747369 

100.00 

90.7 

93.1 

27 

1095009120243 

95.92 

99.9 

100.0 

29 

1095003614744 

100.00 

99.0 

99.8 

30 

1095005227703 

100.00 

98.7 

99.6 
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TABLE  7.1  (Continued) 


Item  no. 

FSN 

Fill  Rates 

Reported^ 

2 

Theoretical 

1  -  8t/Q 

31 

1095105168069 

98.49% 

98.9% 

99.7 

32 

1095004918487 

96.47 

99.7 

99.9 

33 

1095004335657 

77.70 

94.0 

99.3 

34 

1095001111640 

84.73 

79.9 

87.6 

35 

1005008938230 

95.55 

35.3 

94.0 

36 

1005007889718 

99.25 

100.0 

100.0 

37 

1005007379802 

93.55 

97.6 

83.6 

38 

1005007755579 

98.96 

98.6 

99.4 

39 

1005007545293 

100.00 

98.9 

99.7 

40 

1005007545267 

100.00 

96.0 

99.4 

41 

1005007545266 

100.00 

96.9 

99.1 

42 

1005007016793 

100.00 

91.2 

97.2 

43 

1005006999923 

100.00 

39.7 

94.5 

44 

1005006999931 

100.00 

97.9 

99.4 

45 

1 005006999882 

100.00 

91.0 

92.2 

46 

1005000511579 

96.15 

100.0 

100.0 

47 

1005003357318 

100.00 

89.4 

90.5 

48 

1005000178809 

100.00 

100.0 

100.0 

Vrom  fill 

rates  for  April 

1979  -  March  1930 

obtained  from 

LORRA 

2 

From  a  computer  simulation  of  the  current  system 
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Indeed,  there  are  several  methods  of  calculating  theoretical 
fill  ratios.  Brown  (1967,  pp.  91-92)  comments  on  this,  suggesting 


to  us  the  use  of  1-By/Q.  This  alternative  computation  is  also  given 
in  Table  7.1 . 

Methods  and  Analyses  Used 

In  fitting  leadtime  data  to  the  exponential,  the  gamma,  the  normal,  the 
Weibull,  and  the  lognormal,  we  used  an  existing  U.S.  Army  computer  program 
(1971).  In  testing  the  leadtime  data  for  randomness,  we  wrote  our  own 
program.  'Je  similarly  wrote  our  own  computer  programs  in  simulating  leadtime 
demand,  fitting  leadtime  demand  to  the  normal  and  the  Laplace,  doing  ABC 
classification  of  inventory,  and  in  simulating  the  current  and  optimal 
systems.  Furthermore,  we  wrote  a  simulation  program  to  show  that  the  long 
right  tail  of  the  distribution  of  forecast  errors  produced  by  Denny  (1979) 
was  due  to  a  wrong  statistical  procedure. 
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8.  FINDINGS  AND  CONCLUSIONS 

The  major  finding  of  this  study  is  that  leadtime  variability  does  indeed 
influence  safety  stocks,  and  in  order  to  maintain  these  at  appropriate  levels, 
a  one-time  investment  of  about  $150  million  is  necessary.  Also  necessary 
would  be  an  increase  in  annual  operating  expenditure  of  approximately  $50  million. 

Even  if  AFLC  does  not  wish  to  build  up  its  stocks,  it  can  immediately 
halve  its  number  of  backorders  by  explicitly  incorporating  leadtime  variability 
in  its  calculations  and  by  adhering  to  the  Presutti-Trepp  optimal  model.  This 
would  require  a  different  mix  of  items  and,  in  general,  larger  order  quantities. 

We  have  also  found  it  worthwhile  to  use  an  ABC  classification  of  the 
high- intensity  items.  On  the  one  hand,  ABC  classification  is  dynamic  and 
thus  much  superior  to  the  statitic  SMGC  in  current  use  at  AFLC.  On  the  other, 
the  ABC  classification  points  directly  to  the  few  high-intensity  items 
management  should  really  concentrate  its  attention  upon.  In  this  regard 
the  SMGC  seems  anachronistic. 
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APPENDIX  A 

MEMO  OF  UNDERSTANDING 
(MARCH  7,  1980) 
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MEMO  OF  UNDERSTANDING 

SUBJECT:  Hayya  Leadtime  Variability  Study  -  AF  "Contract 

F33615-79-C-5143 

1.  The  results  from  Phase  I  of  AF  Contract  F33615-79-C-5143 ,  Lead- 
time  Variability  in  Inventory  Requirements  Projection,  were 
discussed  in  the  5  Feb  meeting  of  contractor  Dr.  Jack  C.  Hayya,  Capt 
Paul  Gross/AFBRMC  and  representatives  from  LORRA.  During  the  meeting 
it  was  determined  that  Phases  II  &  III  could  not  be  conducted  as 
originally  planned.  The  leadtime  distributions  found  in  Phase  I  were 
for  16  specific  items  during  the  76-79  time  period.  Thus  the  simulation 
model,  IMSSIM,  which  has  71-75  data  cannot  be  used  in  Phase  II. 

2.  Leadtime  data  on  46  more  items  ha^e  been  given  to  Dr.  Hayya  and  he 
has  agreed  to  curve  fit  their  distribution  using  the  same  methods  as 
were  used  on  the  original  16  items.  Also  he  will  determine  and 
analyze  the  convolution  of  the  leadtime  distributions  wich  demand 
distributions  assumed  to  be  Normal,  Laplace,  or  "skewed  Normal".  Dr. 

W.  Steven  Denny's  results  in  "Statistical  Characteristics  of  Forecasting 
Techniques  for  D062  Economic  Order  Quantity  Items"  shows  distribu cions 
for  demands  that  are  near  normal  but  skewed.  The  one  in  Table  IV- 2 
for  OC  with  SMGC=*2  will  be  used  for  this  study. 

3.  A  data  request  for  76-79  demands  for  the  62  items  on  which  Hayya 
has  leadtime  is  being  processed  and  should  be  available  in  March. 

Hayya  will  use  this  to  determine  the  accual  leadtime  demands  for  the 
items  and  test  the  fit  of  the  curves  to  the  p.d.f.s  specified  in  tne 
original  contract.  This  also  will  be  compared  to  the  p.d.f.s  of  the 
convoluted  distributions  discussed  in  paragraph  two  above. 

4.  Dr.  Hayya  has  also  agreed  to  investigate  our  Variable  Safety  Level 
formula.  He  will  conduct  an  analysis  of  the  sensitivity  of  the  safety 
level  to  changes  in  variance  from  demand  variance  to  leadtime  demand 
variance.  Typical  holding  costs  range  from  .15  to  .25.  The  implied 
shortage  factor  or  Lagrangian  multiplier  normally  ranges  from  350  to 
600.  The  Lagrangian  multiplier  is  always  adjusted  so  that  each  ALC 
will  have  a  total  current  safety  level  stock  of  dollar  value  equivalent 
to  2  months  supply  of  all  stock  at  that  ALC.  Any  indication  of  how  the 
mix  of  items  at  an  ALC  will  change  when  using  leadtime  demand  variance 
is  important.  Dr.  Hayya  will  also  make  recommendations  on  possible 
marginal  analysis  techniques  which  will  use  variable  leadtime  demands 
to  minimize  back  orders. 

5.  A  copy  of  "More  Ado  About  Economic  Order  Quantities  (EOQ) "  by 
Victor  J.  Presutti  Jr  and  Richard  C.  Trepp  was  furnished  Dr.  Hayya 
when  he  visited  W-P  AFB  on  5  Feb  1980.  Copies  of  the  EOQ  briefing, 

Vic  Presutti 's  VSL  briefing,  and  "Measurement  and  Implications  of 
Production  Leadtime  Variability",  by  the  US  Army  Inventory  Research 
Office  will  be  furnished  Dr.  Hayya  along  with  a  copy  of  this  memo. 

Demand  data  on  the  62  items  will  be  furnished  as  soon  as  available. 
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EXPLANATION  OF  TABLE  6.1:  OPTIMAL  SYSTEM 
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This  table  can  be  sorted  by  item  no.  or  by  FSN.  The  table  has  24  columns, 
12  on  page  1  and  12  on  page  2.  I  shall  use  the  data  for  item  #48  to  illustrate. 

Col.  1:  Item  no.,  following  the  sequence  in  the  Phase  II  report 

Col.  2:  FSN 

Col.  3:  Daily  Demand,  E(d) 

Col.  4:  Average  leadtime  in  days,  E(L) 

Col.  5:  Variance  of  leadtiine,  V(L) 

Col.  6:  Mean  of  leadtime  Demand 
E(x)  =  E(L)  .  E(D) 

*  (183.25)  (21.70) 

*  3976  ,  for  item  =48 

Col.  7:  Variance  of  leadtime  Demand 

V*(x)  =  E(L)  .  V(D)  +  [E(D)]2  V ( L ) 

=  E(L)  .  E(D)  +  [E(D)]2  V(L)  ,  for  Poisson  demands 
=  E(x)  +  [E(D)]2  V ( L ) 

=  3976  +  (21. 70)2  (535) 

=  255,902  . 

(The  computer  result  of  255,870  is  more  accurate) 

Col.  8:  The  standard  deviation  of  leadtime  demand,  o  (x) 
a  (x)  =^V*(x) 


=  /2558?0 
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Col .  9:  Unit  cost,  ci 

Col.  10:  The  Wilson  Formula 


where  Ai  =  $441.07,  if  c^.  >  $19,500,  $291.82  otherwise 


(See  AFLCR  57-6  (C3),  22  June  1979), 

=  Annual  demand, 

a.  =  holding  cost  factor,  which  is  0.23  in  this  example, 
c.j  =  unit  cost  . 


Hence 


Q 


w 


s 


V  2(441. 07)(650. 96)02) 
(0.23) (3.03) 


*  3113.9  . 


Col.  11:  Optimal  E0Q,  Q  .calculated  according  to  formula  (14)  in  Presutti 
and  Trepp  (1970,  p.  249): 

*  *  (  *2  ) 

Qi  =  0.707  ^  +  |(2AiDi/a.ci)  +  0^/2  > 

★  (2  *2  1 
3  °-707  -i +  {v +  ^ 

*  0.707  (505.84)  +  j(3118.9)2  +  255870/2  |  1/2 


3497  . 


M 


What  is  necessary  in  order  to  use  the  above  formula  is  for  >  2  or  3. 
★  * 

^i  3497 

For  Item  =43,  ttf  ■  50534  *  0.9. 
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* 

Col.  12:  Optimal  safety  factor,  ,  calculated  according  to  formula  IV. 2  in 
Presutti  and  Trepp  (1970,  p.  249): 


Here  we  are  using 

-X  =  660  ,  Z.  =  j  ,  a.  =  0.23  . 

Then 


k1  =  -  0.707  In 


/?( 3497 ) ( 0 . 23 ) ( 3 . 08 ) 


0.5(660)  (0.5)  (505. 84 )(l-e‘^3497^5C5-34^ 


=  -  0.707  In  0.04197 

=  2.242  . 

Col.  13:  Safety  stock  k  j 

kV  «  (2.242))505.84) 

=  1134  ,  for  item  =*48  . 

Col.  14:  Value  of  safety  stock, 

c.kV  =  $3.08  (1134) 

»  $3492  . 


*  -  0.707  1 


_  3503.: 

83,458 
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Col.  15:  Service  level.  This  is 

1  -  \  e^ki*  -  1  -  \  e"*  <2-242> 

=  0.979 
=  0.98  • 


Col.  16:  Weighted  Backorders  at  any  time,  By.  We  may  think  of  this  as  the 

number  of  requisitions  backordered  when  t  1.  If  Z.  =  1,  By  would 
be  the  average  number  of  units  backordered  at  any  time. 

We  use  formula  IV. 2  (Presutti  and  Trepp,  p.  249)  or  formula  (10) 
(Presutti  a. id  Trepp,  p.  246).  From  formula  IV. 2,  we  have 


B 


T 


0.5  Vi 


*2 


For  item  =48, 


0.5 

2 


(})  (255870) 
3497 


(■ 


3497 

505.84 


(0.242) 


=  0.38  . 


The  actual  units  backordered  at  any  time  would  be 


0.38  0.38  . 


0.76  . 
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Col. 


17:  The  fill  ratio  is 


Annual  Demand  -  Annual  backorders 
Annual  Demand 


Annual  backorders  are  given  in  Col.  22.  They  are  calculated  according  ro 
formula  (9)  in  Presutti  and  Trepp  (1970,  p.  246): 


6 


N 


0.5 

/2 


Vi 


★ 


(' 


sz 


k. 


=  —  (-2l-70  x  360)  (505.34)  (0.999943)(o. 041976) 

Ji  "4y/ 

=  16.77  units  per  year  . 


Actually, 


s  *  *  3r  vw 

=  Jl  (0.381(7812)  /  (j) (505. 84) 

=  16.60  (the  difference  is  due  to  rounding)  . 


Hence,  the  fill  ratio  is 
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Col.  24:  Item  ? 
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APPENDIX  C 

EXPLANATION  OF  TABLE  6.2:  AFLC'S  CURRENT  SYSTEM 


•  l 


This  table  can  also  be  sorted  by  item  no.  or  by  FSN.  The  table  has  23 
columns,  11  on  page  1  and  12  on  page  2.  Again  I  use  the  data  for  item  #48  to 
illustrate. 

Col .  1 :  Item  no. 

Col.  2:  FSN 

Col.  3:  Monthly  demand,  E(D) 

Col.  4:  Months  leadtime,  E(L) 

Col.  5:  Leadtime  Demand: 

E(x)  =  E(L)  .  E(D) 

*  (6.11 ) ( 650 . 96 ) 

-  3977,  for  item  no.  48 
Col.  6:  Quarterly  MAD 

Col.  7:  The  std,  deviation  of  leadtime  demand,  according  to  AFLCR  57-6 ( C2 ) , 
29  December  1978.  It  is 

a  *  0.5945  MAD  (0.82375  +  0.42625  L) 


A*.'-, 
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Col.  8:  The  EOQ.  Using  the  EOQ  table  for  WR-ALC  (AFLCR  57-6(C3)  22  June 
1979),  we  have 

Q  =  0.50  ,  since  c^D.  >  $10,000 

=  0. 5(650. 96)(12) 

=  3906  units  . 


Col .  9:  Unit  cost 

Col.  10:  The  safety  factor,  according  to  AFLCR  57-6(C2),  29  December  1978, 
p.  7-2.  It  is 

2/7  Qa.c, 

K  =  -0.707  In  - — -  , 

X(i_)a(i.evrQ/a) 

/R 


where  for  our  example 

R  =  the  average  requisition  size,  4  in  this  case, 
\  =  660, 


and  the  other  terms  as  before.  Hence 


K.  =  -0.707  In 


2/7(3906) (0.23) (3. 08) 

660  ({)  ( 985 . 39 )  (l -e"^ Hot) 


=  -0.707  In 


7826.29 
323 ,983.2 


=  -  0.707  in  0.024156 
*  2.63  . 


The  value  of  K  must  be  between  zero  and  three. 


Col.  11:  Safety  Stock 


ko  a  2.63  (985.39) 

a  2591.58  ,  for  item  #48 


This  value  must  not  exceed  that  for  leadtime  demand. 


Col.  12:  S  Safety  Stock.  This  is 


ciKioi  =  3.08  (2594.27) 
=  $7990.35  • 


Col.  13:  True  K.  The  K-value  of  Col.  10  is  not  based  on  the  true  standard 

★  * 

deviation  of  leadtime  demand,  -  .  The  value  of  j  obtained  from 
Table  6.1,  Col.  8  is  505.84  not  985.39.  Hence  the  true  K  is 


True  K  =  Safety  Stock 


2594.27 

505.34 


5.129  . 


Col.  14:  Actual  3T  weighted  by  essentiality,  Zi .  To  be  consistent,  we  set 


7  «  J_  =  1  . 
1  vR  2 


Then  using  formula  IV. 2  in  Presutti  and  Trepp  (1970,  p.  249),  we  have 

•,•¥*£  ('--4) 

=  ¥(|hW  (, .  e-.^)) 


0.006 
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The  actual  number  of  backorders  expected  at  any  time  is  actually 


8t/Z.  =  0.006/0.5  =  0.012  . 

Col.  15:  Actual  Service  Level 

* 

1  -  1  e”*  ki 
=  1  -  (0.5)  (0.00070769) 

=  1.00  • 


Col.  16:  $  Inventory  on  hand.  This  is 

ci  (Safety  stock  +  |  +  3T/Z.) 
»  3.08  (2594.27  +^^+0.012) 
=  514,005. 

Col.  17:  Fill  Ratio.  This  is 


Annual  Demand  -  Annual  Number  of  backorders 
Annual  Demand 


We  obtain  the  annual  number  of  backorders  using  formula  (9)  in  Presutti- 
Trepp  (1970,  p.  246): 


8 


N 


0.5 

/I 


(■ 


=  02)  (650^96)  (50S. 84)  (,3  99993)  (q. 00070769) 

/2" 


0.25  . 


h*  ' 


You  see  this  value  in  Column  21. 


Hence, 


ri'.i  „_t4-  _  12(650.96)  -  0.25 
Fill  Ratio  =  — i -rm. - 


=  1.00  - 


Col.  18:  Annual  Ordering  Cost 


441 .07(650.96) (12! 

im 


=  $882.09  =  $882.14 


It  is  strange  that  the  ordering  cost  for  so  many  items  is  $882.14  or 
$583.64.  Apparently  an  annual  ordering  cost  of  $382.14  corresponds  to  an 
ordering  cost  of  S441.07,  whereas  an  annual  ordering  cost  of  $583.74  corresponds 
to  an  ordering  cost  of  $291.32.  Thus 

$HtTW'  (5882.14)  =  $583.64  . 

Consider  item  =3  where  the  dollar  value  of  annual  demands  is  373.92  X 
(12)($2.74)  =  $12,294.49  <  $19,500.  Hence,  the  ordering  cost,  A  *  $291.82; 
and  the  annual  ordering  cost  would  be 


A0  .  (291 .82 ) (373. 92) (1 2 1 
T  "  2245 


»  $583.52  . 

It  seems  that  the  E0Q  table  for  WR-AIC  (AFLCR  57-6(C3),  22  June  1979, 
p.  7-14-1)  is  designed  to  produce  precisely  these  annual  order  costs. 
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Col.  19:  Annual  Holding  Cost  based  on  the  Inventory  Position.  This  is 
a^ci  (y  +  Safety  Stock  +  |  ) 

=  $0.23(3.03) (3976  +  2594.27  +  ) 

=  $6037.88  . 

Col.  20:  Annual  operating  cost.  This  is  the  sum  of  the  annual  ordering  cost 
and  the  annual  holding  cost. 

Col.  21:  Annual  number  of  units  backordered.  This  has  already  been  calculated 
under  Col.  17. 

Col.  22:  FSN 


Col .  23:  Item  no. 
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